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ABSTRACT 



Purpose: The aim of this study was to provide sex-matched three-dimensional (3D) statistical shape models of the 
mandible, which would provide cephalometric parameters for 3D treatment planning and cephalometric 
measurements in orthognathic surgery. 

Materials and Methods: The subjects used to create the 3D shape models of the mandible included 23 males and 
23 females. The mandibles were segmented semi-automatically from 3D facial CT images. Each individual 
mandible shape was reconstructed as a 3D surface model, which was parameterized to establish correspondence 
between different individual surfaces. The principal component analysis (PC A) applied to all mandible shapes 
produced a mean model and characteristic models of variation. The cephalometric parameters were measured 
directly from the mean models to evaluate the 3D shape models. The means of the measured parameters were 
compared with those from other conventional studies. The male and female 3D statistical mean models were 
developed from 23 individual mandibles, respectively. 

Results: The male and female characteristic shapes of variation produced by PC A showed a large variability 
included in the individual mandibles. The cephalometric measurements from the developed models were very close 
to those from some conventional studies. 

Conclusion: We described the construction of 3D mandibular shape models and presented the application of the 3D 
mandibular template in cephalometric measurements. Optimal reference models determined from variations pro- 
duced by PCA could be used for craniofacial patients with various types of skeletal shape. (Imaging Sci Dent 2012; 
42 : 175-82) 
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Introduction 

In modern surgical treatments, three-dimensional (3D) 
computed tomography (CT) plays a fundamental role in 
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estimating the final treatment outcome and selecting the 
appropriate surgical intervention. 1 ' 2 Three-dimensional CT 
can provide precise and detailed information for the diag- 
nosis of craniofacial structural problems and for the preop- 
erative simulation of operations such as treatment of cran- 
iofacial anomalies and bone defects, and trauma surgery. 3 " 6 
Patients with distinct craniofacial deformities or missing 
bony structures require surgical reconstruction. In planning 
craniofacial surgical interventions, the symmetrical appear- 
ance in the frontal view and balanced facial profile in the 
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lateral view as well as the functional aspects of the patient 
are important considerations. When those cannot be achi- 
eved by a mirror image of the unaffected contralateral side, 
a three-dimensional reference model of craniofacial bone 
is necessary to compare the pathologic conditions with a 
normal anatomy and finally to construct a harmonized 
facial appearance. Zachow et al 7 developed a statistical 
3D mandibular model from 1 1 different persons and used 
this model as a template for the planning of surgical recon- 
struction in cases of mandible deformities. Gateno et al 8 
reported that the computerized composite skull models 
were created and used for computer-aided surgical simula- 
tion in treatment of complex craniofacial deformities. 
There-dimensional medical modeling methods based on 
CT data have been widely used in planning orthognathic 

9,10 

surgeries. 

3D cephalometric analysis is essential for the computer- 
assisted planning of craniofacial surgical procedures, par- 
ticularly in orthognathic surgery. The three-dimensional 
cephalometric parameters measured from 3D reconstruct- 
ed images have been used in the diagnosis of facial asym- 
metry, 1112 in the simulation workbench for orthognathic 
surgery, 13 and in the evaluation of postoperative changes 
of mandibular anatomy and position after orthognathic 
surgery. 14 In those studies, the cephalometric parameters 
were averaged measurements from the individual refor- 
matted 3D images, which only provided the specified 
parameters. Recent orthognathic surgical techniques and 
treatments require more advanced and complex informa- 
tion from diagnostic images. 15 Three-dimensional cephalo- 
metric analysis providing optimal references for patients 
of various skeletal types would be useful for clinical appli- 
cation of computer-aided techniques in craniofacial surgery. 
The aim of this study was to provide sex-matched statisti- 
cal shape models of the mandible, which would provide 
cephalometric parameters for 3D treatment planning and 
cephalometric measurements in orthognathic surgery. 
Three-dimensional statistical shape models of the mandible 
were developed using principal component analysis (PC A). 
The 3D cephalometric parameters were directly measured 
from the mean models and compared with those from other 
conventional cephalometric analyses. 

Materials and Methods 

The subjects consisted of 23 males (mean age, 24.7 years) 
and 23 females (mean age, 26.0 years) from our dental 
hospital. All subjects had Angle class I molar relationship 




Fig. 1 . Decomposition of 3D mandible shape for the construction 
of correspondence maps. 

without mandibular asymmetry recognizable to the naked 
eye. The CT images were obtained using a Somatom Sen- 
sation 10 (Siemens, Erlangen, Germany) with the same 
imaging condition (slice thickness: 0.75 mm, slice interval: 
0.5 mm, 120 kVp, 100 mAs). 

The development of the 3D statistical shape models of 
the mandible consisted of four steps: segmentation of the 
mandibles, construction of correspondence maps, align- 
ment, and principal component analysis. In the segmenta- 
tion step, the mandible was semi-automatically segmented 
from the 3D facial CT images by a dentist with clinical 
experience in the Department of Oral and Maxillofacial 
Radiology. The segmentation procedure included 3D 
thresholding, 2D region growing, 3D region growing, mor- 
phological operations, and contour tracking. After segmen- 
tation, the 3D individual mandible shape was reconstructed 
as a surface model with triangulated meshes. Smoothing 
and simplification of the meshes were applied to the shape 
model. 

For construction of correspondence maps, each 3D 
mandible surface was decomposed into 8 corresponding 
patches in a symmetric manner (Fig. 1). Then, the mandible 
was split in half and subdivided through the lower central 
incisors. The condylar head, ramus, and mandibular body 
were separated on each side. The tooth region was excluded 
to separate variation in individual dentition. Then, each 
corresponding patch was mapped to a common disk using 
homeomorphic mapping under minimal geometric distor- 
tion. That is, it was parameterized to establish correspon- 
dence between different individual surfaces, which yielded 
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Table 1 . Landmarks used for 3D cephalometric measurements in different viewing positions 





Landmarks 


Description 


Anterior view 


Menton (Me_ ant ) 


The most inferior point on the symphyseal outline in anterior view 


Posterior view 


Condylion posterioris (Cd post post ) 
Gonion posterioris (Go post post ) 


The most posterior point of the condyle head in posterior view 

The most posterior point of curvature along the angle of mandible in posterior view 


Superior view 


Condylion superioris (Cd sup . sup ) 


The most superior point of the condyle head in superior view 


Inferior view 


Gonion inferioris (Go inf inf ) 

/A.lllCgUlllUll 11UIC11 v^/^.g inf/ 

Menton (Me inf ) 


The most inferior point of curvature along the angle of mandible in inferior view 

TVif^ mnst inffrinr 'nnint nf anterior Vinrrlfr nf tVif ancrlf^ nf manrlihlf in inf print* vipw 

111C lllUaL 1111C11U1 JJUllll Ul CllllCllUl UU1U.C1 Ul LUC dllglC Ul llldllU-lUlC 111 1111C11U1 V1CW 

The most inferior point on the symphyseal outline in inferior view 


Lateral view 


Condylion posterioris (Cd post _ lat ) 
Gonion posterioris (Go post _ lat ) 
Gonion inferioris (Go inf lat ) 
Antegonion(Ag i at ) 
Sigmoid notch (S) 


The most posterior point of the condyle head in lateral view 

The most posterior point of curvature along the angle of mandible in lateral view 

The most inferior point of curvature along the angle of mandible in lateral view 

The most superior point of antegonial notch of mandible in lateral view 

The most inferior point of sigmoid notch 




Goi 

Go 2 

Go 3 

Gn 

Cd 

Cp 


The most posterior point of posterior border of ramus 

Midpoint of posterior border of ramus 

The most inferior point posterior border of ramus 

The most anterioinferior point on the contour of the symphysis 

The most superior point on the head of the condyle 

The most upper and posterior aspect of condyle 



the correspondence maps. The subdivided surfaces pre- 
served all the topological properties of a disk while estab- 
lishing correspondence maps. 

In an alignment step, first, the shape with the highest 
number of mesh points was selected as the reference of the 
mandible shapes. Registration of a shape as a reference for 
all of the other shapes based on a rigid-body transforma- 
tion was achieved by minimizing the sum of the squared 
distances between the corresponding points. After these 
processes, principal component analysis (PCA) was simul- 
taneously applied to all the patches of the mandible by 
achieving continuity across patches using the common 
boundaries used for decomposition. All mandible shapes 
were represented as vectors in the 3N dimension (N, the 
number of reference mesh points). First, the mean shape 
was computed by averaging the vectors of all of the differ- 
ent shapes. Then, the characteristic models of variation 
were generated by a linear combination of the eigenmodes 
produced by PCA. 7 ' 16 

The principal component analysis transformed a number 
of possibly correlated variables into a smaller number of 
uncorrelated variables called principal components. The 
PCA produced the eigenvalues (Ai> A 2 > ... > A^O) and 
the corresponding eigenvectors (p h 2, n). Only 
the few first modes carried the most important informa- 
tion; therefore, each original shape was reconstructed using 
only some principal components. As a result, it was possi- 
ble to represent as much variance contained in the train- 



ing set as possible by as few parameters as possible. 

First, the fundamental cephalometric parameters were 
measured from the mean models to evaluate the 3D shape 
models that had been developed. Eighteen landmarks 
according to the different viewing positions are explained 
in Table 1 and illustrated in detail in Figure 2. To compare 
the measurements with those from various viewing posi- 
tions by other studies, multiple landmarks were selected 
for a single anatomical region. Using 3D measurement 
functions, two dentists directly measured 23 cephalometric 
parameters from male and female models, respectively. In 
comparable studies, the measurements were performed 
from the individual reformatted 3D images of the CT scan 
data. Therefore, the cephalometric parameters from these 
reports were averaged values of individual measurements. 
The nineteen parameters were compared with the results 
of Ahn et al, 17 which were measured from the mandibles 
of 60 normal occlusion individuals (30 males and 30 fe- 
males) for the diagnosis of facial asymmetry. The three 
parameters were compared with the results of Park et al, 18 
which were measured for the analysis of craniofacial mor- 
phology (16 males and 14 females). The nineteen parame- 
ters used for comparing with the results of Ahn et al 17 are 
as follows: 

(1) Ramus length: Cd sup sup - Go in f Jat , Cd sup _ sup - Go inf _ inf , 
Cd sup _ sup — Go pos t_iat? Cd sup _ sup — Go pos t_ po sb Cd sup _ sup — 

Agiat, Cd sup _ sup — Agi n f, 

(2) Mandibular body length: Go inf _i at - Me_ ant , Go inf _i at - 
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Fig. 2. Established landmarks for 3D cephalometric measurements 
in different viewing positions. 

Me_i n f, Goinfjnf— Me_ ant , Goinfjnf— Me_i n f, Go post j at — 
Me_ ant , Gop OSt j at — Me_i n f, Gop OSt _p OSt — Me_ ant , Gop OSt _p OSt — 
Me_ inf , 

(3) Condylar neck length: Cd sup _ sup — S, 

(4) Gonial angle: ZCd post _i at -Go post _i at -Me_ ant , ZCd posUat - 
GOpostjat — Me_i n f, Z_Cdp OSt _p OSt - Gop OSt _p OSt - Me_ ant , 
^-Cdp OSt _ post — Go post _p OSt - Me_i n f . 

The three parameters compared with the results of Park 
et al 18 are as follows: 

(1) Ramus length: Cp-Goi+Goi-Go 2 

(2) Mandibular body length: Me-Go 3 +Go 3 -Go 2 

(3) Gonial angle: Cp — Go 2 — Me. 

The one parameter compared with the results of Kim et 
al 19 is as follows: 

(1) Mandibular body length: Cd-Gn. 

In the comparison with other studies, the statistical infor- 
mation we were able to obtain from the other studies only 
included the means of the measurements. Therefore, the 
most reasonable analysis was a simple comparison using 
the differences between the means. The means of the mea- 




Fig. 3. Male and female 3D statistical mean models (left: male, 
right: female). 

surements by two observers differed from those of other 
studies. 

Results 

Male and female statistical mean models were created 
from each of the 23 individual mandibles. The male and 
female 3D statistical mean models are shown in Figure 3. 
The male and female characteristic shapes of variation pro- 
duced by PC A are also shown in Figures 4 and 5, respec- 
tively. Figure 4 shows the reconstructed shapes by varying 
the first three modes of variation for the male mandibles. 
In the first row, the mode corresponding to the largest vari- 
ance (Ai) is varied from -3i/A~i to 3\[X\ , and in the second 
and the third row, the modes correspond to the second and 
the third mode, respectively. Figure 5 shows the recon- 
structed shapes by varying the first three modes of varia- 
tion for female mandibles. The statistical shape models 
show the large variability included in the individual mandi- 
bles. 

The cephalometric parameters were calculated as the 
average of measurements determined by two dentists using 
the mean models. The results from the male and female 
mean models are shown in Table 2. All the cephalometric 
parameters measured in the male were greater than those 
in the female except ZCd post _i at -Gop OSt _i at -Me_ in f and 
^Cdp OSt _p OSt -Go post _p OSt -Me_ in f of the gonial angles. The 
differences between measurements from the male mean 
model and the results of Ahn et al 17 were 0.7-4.3 mm at the 
ramus length, 0.3-4.1 mm at the mandibular body length, 
and 0° -1.7° at the gonial angle. In the female model, the 
differences were 1.7-2.7 mm at the ramus length, 0.2-2.4 
mm at the mandibular body length, 0.2 mm at the condylar 
neck length, and 1.9° -3.4° at the gonial angle (Figs. 6-8). 
The differences between measurements from the male 
model and the results of Park et al 18 were 3.1 mm at the 
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b=-3/I b=-2/I b=-JX b=0 b=JJ b=2/X b=3/X 

Fig. 4. Characteristic models of variation produced by PCA for male mandibles. The mode corresponding to the largest variance (Ai) is 
varied from to -3{X\ to 3{Xi (the first row), and the modes corresponding to the second mode (the second row) and the third mode (the 
third row). 



b=-3jX b=-2/I b=-JX b=0 b=JI b=2pL b=3jl 

Fig. 5. Characteristic models of variation produced by PCA for female mandibles. The mode corresponding to the largest variance (Ai) is 
varied from to -3{X\ to 3{X\ (the first row), and the modes corresponding to the second mode (the second row) and the third mode (the 
third row). 



ramus length, 2.0 mm at the mandibular body length, and 
0.7° at the gonial angle. In the female model, the difference 
values were 7.5 mm at the ramus length, 0.5 mm at the 
mandibular body length, and 6.4° at the gonial angle (Figs. 
6-8). The difference between measurements from the male 
model and the results of Kim et al 19 was 0.7 mm at the 
mandibular length. In the female model, the difference was 



0.2 mm (Fig. 7). The results from the developed models 
were very close to those from other conventional studies. 

Discussion 

In cases of craniofacial malformations and acquired de- 
fects, the defected and malformed structures are recon- 
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structed based on normal anatomic structures. When the 
affected structures are not complex, the affected half of 

Table 2. Cephalometric measurement results from the male and 
female mean models 

Measurements parameters Male model Female model 



the structure can be replaced by a mirror image from the 
unaffected side. 12 However, in complex cases, the recon- 
struction cannot be guided by mirror images of the unaf- 
fected side. A surgeon must compare the pathologic situa- 
tion with a mental image of normal anatomy. 20 In those 
cases, a 3D optimal shape model can be used as a reference 
template for the anatomical region. The shape models of 
variation that we developed would be able to serve as 3D 
templates for the reconstruction of missing or malformed 
mandible structures. 

The 3D cephalometric analysis is essential for computer- 
assisted planning and interventions in orthognathic surgery. 
Three-dimensional CT measurements have been applied 
in treatment of craniofacial malformations and acquired 
defects, as a planning tool for the skull base reconstructive 
surgery, 21 for surgical planning of head and neck cancer, 22 
and for 3D landmark measurement in craniofacial surgery 
planning. 23 The collected CT data has been sorted and 
classified according to sex and age to serve as sex- and age- 
dependent norm data, which can be used as a basis for vir- 
tual patient-specific operation planning and simulation. 20 
While the cephalometric parameters are averaged measure- 
ments from individual reformatted 3D images in those 
studies, the 3D cephalometric parameters can be measur- 
ed from the mean models directly in this study. 

The cephalometric parameters measured from the mean 
models were compared with those of other previous studies 
to evaluate the 3D mandible models that we developed. 
The differences in the male model were 0.3-4.3 mm in 
length and 0° -1.7° in angle, and the differences in the fe- 



Ramus length (mm) 

Cd SU p. SU p GOi n f_i a t 

Cd SU p_ SU p Goj n f_i n f 

Cd SU p. SU p GOp OS t_iat 

Cd SU p_ SU p Gop OS t_p OS t 
Gd SU p_ SU p Agj at 
Gd SU p_ SU p Agj n f 
Cp-Goi+Goi-Go 2 

Mandibular body length (mm) 
Go in f_iat-Me_ ant 
Go in f_ lat -Me_ in f 
Go in f_ in f-Me_ ant 

GOi nf j n f-Me_i nf 

Go pos t_ lat -Me 

ant 

Go p ost_lat — Me_i n f 

GOp OS t_post — Me_ ant 
Go post _p OSt -Me_i nf 
Me - Go 3 +Go 3 - Go 2 
Cd-Gn 

Condylar neck length (mm) 

Cd SU p_ SU p — S 

Gonial angle (°) 

^-Cdpostjat — Gop OS tjat — Me_ an t 

^-Cdpostjat — GOp OS t_iat — Me_i n f 

-^-Cdp OS t_post — Gop OS t_p OS t — Me_ ant 

^-Cdpostpost — GOp OS t_post — Me_i n f 

ZCp-Go 2 -Me 



69 


6 + 4 


18 


65 


4 + 5 


12 


70 


4 + 4 


26 


66 


3 + 4 


39 


55 


0 + 2 


77 


54 


1 + 3 


29 


54 


4 + 3 


42 


52 


9 + 2 


10 


74 


1 + 3 


86 


69 


6 + 3 


67 


74 


6 + 5 


28 


70 


2 + 5 


14 


63 


4 + 5 


72 


58 


5 + 3 


70 


83 


3±4 


91 


81 


5 + 3 


72 


82 


2 + 3 


69 


79 


1 + 3 


17 


82 


6 + 4 


94 


79 


2 + 2 


09 


81 


5 + 3 


08 


78 


2 + 2 


40 


96 


1+4 


82 


94 


6 + 4 


40 


94 


4±4 


20 


92 


0 + 3 


06 


96 


4 + 3 


18 


94 


3 + 3 


83 


93 


4 + 4 


40 


92 


6 + 3 


51 


93 


7 + 4 


73 


89 


7 + 3 


83 


129 


1 + 5 


50 


121 


4 + 4 


50 


27 


5 + 1 


24 


26 


2+1 


30 


121 


3±4 


46 


120 


8 + 3 


54 


121 


7 + 5 


43 


122 


9 + 3 


25 


122 


0 + 5 


39 


121 


5 + 4 


28 


120 


8 + 4 


74 


122 


5 + 4 


39 


119 


5 + 4 


65 


117 


9 + 2 


88 




Cd SU p_ SU p Cd SU p_ SU p Cd SU p_ SU p Cd SU p_ SU p Cd SU p_ SU p Cd SU p_ SU p Cp G01+G01- 

GOi n f_i at GOi n f_i n f Go p0 st_lat Go post _p OSt Agi at Agj n f GO2 



Ramus lengths 
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Mandibular body lengths 




^Cdp 0S t j at - ZCdp 0S t j at - ZCdp 0S t post - ZCdp 0St p 0St - Cp-Go 2 -Me 

GOpostjat - Me_ ant Go post jat ~ Me _ inf Go post _ pos t - Me_ ant GOp OS t_p OS t - Me _ M 

Gonial angles 

Fig. 8. Comparison of gonial angles from the developed models 
with Ahn et al 17 and Park et al. 18 



male model were 0.2-2.7 mm in length and 1.9° -3.4° in 
angle. The data measured from the models were very close 
to those from other conventional studies. The differences 
in the sample composition and measurement methods 
might explain the differences. In order to compare with 
other studies that only provided a specified number of 
parameters, 23 parameters were measured from the mean 
models that we developed. Since the developed shape 
models can be observed and cephalometric parameters can 
be conveniently measured from any direction, it is possible 
to provide more diversity in the 3D cephalometric mea- 
surements. Therefore, other cephalometric measurements 
can be performed to apply diverse 3D cephalometric para- 
meters in maxillofacial surgery planning and orthodontic 
treatments. In this study, the cephalometric parameters 
were only measured on the right side of the mandible. The 



measurements on both sides of the mandible models can 
be performed to evaluate the mandibular asymmetry in 
further studies. The mean models developed in this study 
are not sufficient to be called the standard model for the 
time being. By including more normal mandible CT data 
in developing the mean models, the models will more 
accurately approximate the Korean normal 3D mandible. 
The measured parameters will also approach the standard 
value of the Korean norm. 

Optimal reference models for patients can be determin- 
ed from the variation models that best match the measure- 
ments from the healthy part of the mandible. New cephalo- 
metric parameters can be measured based on the evaluation 
of skeletal landmarks on the models of variation produced 
by PCA. These parameters will be useful in the diagnosis, 
surgical planning, and follow-up after surgery of cranio- 
facial patients with various types of skeletal shape. 

When a surgeon plans computer-aided craniofacial sur- 
gery, it is generally necessary to segment and reconstruct 
the individual mandible of patients from 3D CT images. 
The streak artifact caused by the metal restorations in teeth 
is a frequently encountered problem with 3D CT maxillo- 
facial imaging. The editing process to delete the metal 
artifacts from each axial slice of the image is very time- 
consuming and labor-intensive. This statistical shape mod- 
el can be used for automatically segmenting the 3D anatom- 
ical structures. 16 The model-based segmentation using a 
3D mean shape model can provide more efficient segmen- 
tation of the mandible in dental clinics. For developing the 
3D shape model, we used the method of patch decompo- 
sition and parameterization to establish correspondence 
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between the different individual surfaces and principal 
component analysis to yield mean and variation shapes. 
The same method and analysis could also be applied to 
develop other statistical 3D mean models such as for the 
maxilla, zygomatic bone, and teeth. 

We described the construction of 3D mandibular models 
and presented the application of the 3D mandibular tem- 
plate in 3D cephalometric measurements. Optimal refer- 
ence models determined by the variations produced by 
PCA can be used for craniofacial patients with various 
types of skeletal shape. Due to the recent widespread in- 
crease in the use of cone beam CT, the acquisition of 3D 
craniofacial data is no longer difficult. As a result, the de- 
mand for standardized 3D models will also increase in den- 
tal clinics. 
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